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The probability distribution of the current in the asymmetric simple exclusion process is expected 
to undergo a phase transition in the regime of weak asymmetry of the jumping rates. This transition 
was first predicted by Bodineau and Derrida using a linear stability analysis of the hydrodynamical 
limit of the process and further arguments have been given by Mallick and Prolhac. However it has 
been impossible so far to study what happens after the transition. The present paper presents an 
analysis of the large deviation function of the current on both sides of the transition from a Bethe 
Ansatz approach of the weak asymmetry regime of the exclusion process. 



I. INTRODUCTION 



The simple exclusion process is one of the toy-models of out-of-equilibrium statistical physics and was introduced 
first around forty years ago pQ. Many aspects of equilibrium statistical physics of particle systems are now well 
understood, however only a few approaches are known for out-of-equilibrium systems. Such systems usually have a 
non-trivial irreversible dynamics that prevents the stationary regime to be Boltzmann-Gibbs-like. Moreover, they are 
also characterized by the existence of non-zero currents (energy, mass, etc.) going through the system. 

The asymmetric simple exclusion process (ASEP), as defined below, is such a model of mass transport along a 
one-dimensional periodic track of length L: each site is occupied by at most one particle and each particle jumps to 
the right (resp. left) with a rate p (resp. q) if the target site is empty. The exclusion rule induces an interaction 
between neighbouring particles. If p ^ q, a current of particles flows through the system. 

Close to equilibrium, for example if the difference p — q cx 1/L is small on a large lattice, the system can reach local 
equilibrium and it can be described at the macroscopic scale by a hydrodynamical theory [2-6 . This theory allows in 
principle to compute the typical macroscopic density profile that corresponds to the observation of a given current 
as well as the large deviation of this current from a variational principle. Finding the global optimal profile in this 
framework remains however a difficult task [7] [8] : the flat density profile is the optimal profile for currents close to the 
mean current but a linear stability analysis shows that this optimal profile becomes unstable and time-dependent profile 
for smaller currents. Although the scenario of a traffic jam occurring at the transition seems clear, the characterization 
of the transition and of the jammed profile was still unknown. 

A second approach to the asymmetric exclusion process comes from the integrability of the microscopic model and 




Figure 1: Asymmetric exclusion process on a ring geometry of L = 16 sites with N = 6 particles 
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the relation between its Markov matrix and the Haniiltonian of the XXZ quantum spin chain. Formally, the Markov 
matrix modified so that the current can be counted can be diagonalized with the Bethe Ansatz and the eigenvalues 
parametrized by complex numbers satisfying coupled non-linear equations. This approach has already lead to many 
results about the stationary state of exclusion process and the slowest relaxation times on the ring [9"HT4] and more 
recently with couplings to two reservoirs |15fH5] . It has been also possible to study the fluctuations of the current 
both numerically [19 , 20 and theoretically through the systematic perturbative computation of the cumulants of the 
current [2T] l2"3~H2l)] . The present paper shows how to go beyond these results and describe the jammed phase as well as 
the transition itself from the point of view of the Bethe Ansatz. 

The paper is organized as follows: section [II] introduces the model, the notations and some previous results. In 
section |III| the Bethe equations are rewritten in a suitable way to take the simultaneous limit of large size and weak 
asymmetry. The Gaussian part of the probability distribution of the current is recovered in section |IV A| from the 
Bethe Ansatz, the transition in section [lVB| and the jammed phase in |IV C[ followed by a more general discussion of 
the type of phase transition that is described. 



II. GENERAL PROPERTIES 



A. Definitions 



The asymmetric exclusion process describes N particles evolving on a one-dimensional lattice of length L with 
periodic boundary conditions (see figure [lj . Each site i is occupied by at most 1 particle and is described by an 
occupation number 7$ £ {0, 1}. The positions of the particles are noted Xj for 1 < j < N so that Tj — 1 if and only if 
i = Xj for some j. 

During a infinitesimal time dt, a particle on a site i tries to jump with rate p to the site i + 1 and with rate q to the 
site i — 1: if the target site is empty, the jump occurs, whereas it is cancelled if it is occupied. In the two-site basis 
(00, 01, 10, 11), the Markov matrix that describes the jumps through the corresponding edge is 

/0 0\ 
\0 0/ 

The Markov stochastic matrix W for the whole lattice is then given by 

L 

W = £>.•,<+!> (2) 
1=1 

where the site L + 1 is identified to the site 1 and WJj j+i is the matrix w above acting on the pair of sites (i, i + 1). 
This dynamics keeps the number of particles constant and one can label each sector of W by the density po — N / L. 

A slight modification of the matrix W gives the Legendre transform of the large deviation function of the current. 
The current Q accumulated between and t is defined as the difference between the total numbers of jumps of particles 
to the right and of jumps to the left. The probability P t (C, Q) of observing the configuration C at time t and a total 
current Q between and t satisfies 

8 t P t (C, Q)=J2 W ccP t (C, Q ~ qc>,c) - W ^c F t (C, Q) , (3) 

where Wqc is the matrix element of VV which gives the transition rate from C to C and qc e € {—1, 1} is the increment 
of the current corresponding to the change from C to C. The generating function of the current Pt (C, s) defined as 

P t (C,s) = ]Te s< 3p t (C,Q) (4) 



satisfies a modified version of (§: 



d t P t (C,s) = (e S9c ' c 'Wcc<)P t (C',s) - I Wcc J 



P t (C,s). (5) 
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This equation is linear and its linear operator yy( s ) is given by 



W« = X>m + i< ( 6 ) 

/0 0\ 

-g pe s ^ 



(?e~ s -p 
Vo 0/ 



The matrix W^ s ^ is not stochastic any longer but its eigenvalue Ai(s) with the maximal real part is real and 
non-degenerate (Perron- Frobenius' theorem) in each sector N. This eigenvalue gives the long time behaviour 

Pi (C,s) oa (f>[ s \c)e Xl ^ t where cf)±(C) is the associated eigenvector. This exponential growth corresponds to a large 
deviation behaviour of the distribution of the current: 

where f(J) and Ai(s) are related through the Legendre transform 

A 1 ( S )=su P (/(j) + SJ ). (9) 

3 ' 

The expected phase transition described in the introduction |S] occurs in the function f(j) (resp. Ai(s)) as the 
parameter j (resp. s) varies. The function /(j) is negative and is zero when j is equal to the mean current. The 
successive derivatives of Ai(s) at s = give the cumulants of the current (Q") c already given for some cases in 
[H[21[25]. 

Once normalized, the eigenvector (C) describes the stationary measure of the configurations conditioned on the 

past current. The matrix has left eigenvectors {C) associated to the same eigenvalues as the right eigenvectors 

4>\ S \C). For s = 0, the stochasticity of W implies that the eigenvector associated to Ai(0) = is given by ipf\C) = 1 

up to a normalization constant. As explained for example in |18| . the knowledge of ip[ s \c) is necessary to understand 
the distribution of the configurations conditioned to the production of a given current j at larger times. The probability 
of observing a configuration C at time t — aT for some a given in ]0, 1[ knowing that the current at the final time T is 

equal to Qt = jT converges as T — > oo to V^J nd (C) given by: 

^oU C ) = ^i i) (C)4 Sj) (C) (10) 

where Sj = —f'(j) and Z is a normalization constant (see |18|). 

Symmetries. The asymmetric exclusion process is invariant under the particle-hole transformation. For a system of 
N particles, the L — N holes jump to the left with rate p, to the right with rate q and satisfy the exclusion rules. The 
system is also translation-invariant. A less trivial symmetry is the Gallavotti-Cohen symmetry [27] that relates the 
probabilities of time-reversed histories of the system. One can check on the definition of that it satisfies: 

W (s)T = W (s '\ with s + s' + ln(p/g) = (11) 

where M T denotes the transposition of the matrix M. It follows that the first eigenvalue Ai(s), which is non-degenerate, 
satisfies: 

X 1 ( S ) = X 1 (-\n(p/q)- S ). (12) 

This identity in terms of the large deviation function f(j) can be written f(j) — f(—j) = ln(p/q)j. 

Scalings. The following sections give the largest eigenvalue Ai(s) for any s in the limit of weak asymmetry defined 
as: 

(13a) 
(13b) 
(13c) 



p = 


1 


V 

ho 

2L 




q = 


1 

2 " 


V 

ho 

2L 




s = 


7 

L ' 


fo^ 1 ) 
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This scaling corresponds to the scaling already used in many different works (see [7J H3] for examples) but gives very 
different results from the scaling p — q oc 1/y/L, which is often also called the weak asymmetry regime and is related 
to the KPZ universality class [2"5H3"U] . 

Remark that this scaling is not usual in spin chain models such as the XXZ spin chain, for which the coupling 
constants are not assumed to depend on the size of the system. However, in the present case, the scaling is natural 
from the point of view of the discretization of a macroscopic continuous diffusion process. 



B. Hydrodynamical description 



In the weak asymmetry regime ( 13 1 , the system can be described on the macroscopic scale by the hydrodynamic 
limit studied in [21|7j|8]. A macroscopic configuration is characterized by its density profile p(x,t) at time t where 
x G [0, 1] is the rescaled position. The large deviation function f(ja) corresponding to the observation of a current jo 
over a large time T is given in this approach by a variational principle 



(14) 



where the local current j(x,r) is an auxiliary variable which must satisfy both the global current constraint 



and the local conservation of the mass 



1 f T 

Jim ™ / j{x,T)dT = j (15) 



d T p{x, t) + d x j(x, t) = 0. (16) 



The position- and time-dependent density profile must also satisfy the conservation law p(x, r)dx — p. The functional 
I\0,T] (jy p) is given for the weakly asymmetric simple exclusion process by: 



T f 1 (j{x,T) + D(p(x,T))d x p(x,T)-va(p(x,T))) 



2 

%T]0',P) = / dr I dx± — — (17) 

Jo Jo 2 <r(p{x,V) 

D(p) = \, (18) 
a(p) = p(l-p). (19) 

The two functions D and a given in this equation are particular to the WASEP but the expression of I[o,t] {j, p) is 
more general and can describe other non-equilibrium diffusive systems [21 [5] . 

Minimizing (14 1 is not easy since the optimal profile may become time-dependent. For j close to the mean current, 
the flat solution p(x,r) = po is found to be the optimal profile [7 . In this case, a simple computation shows that the 
distribution of the current is nearly Gaussian and thus one has Lf(jn) — > —(jo — is<t(pq)) 2 / (2a(po)) and 

L Al ( S = 7 /i) ^ Af uss ( 7 ) = ^1 ( 7 2 + 2vi) . (20) 

A phase transition is however expected to occur when j = j c (and 7 = j c for Ai(7)) since it has been shown in [7| 
that after j c the flat profile is linearly unstable for the evolution induced by X[ 0j T](jiP)- The critical parameter 7 C is 
obtained by solving 

Af uss ( 7 ) = -n 2 /2. (21) 

For Af auss (7) < — 7r 2 , the characterization of the optimal profile is still missing although candidates with a travelling 
wave shape p(x,r) — g(x — vt) have been studied 0[S] and found to produce a lower value of T[n t T]{j> p) than the 
flat profile. The minimization of (17 1 over the moving profiles g on the ring and the velocities v leads to a profile g 
determined only implicitly from a set of elliptic integrals. The details of the equations that are obtained are presented 
in appendix |A") 
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C. Integrability of the microscopic dynamics 



From an algebraic point of view, the matrix can be mapped onto the Hamiltonian of an XXZ spin chain with 
periodic twisted non-hermitian boundary conditions, which is known to be integrable (see [9l412( 13 1 j for earlier results 
on ASEP-like models based on Bethe Ansatz methods). The only result from integrability used in this section is that 
the matrix can be diagonalized with the coordinate Bethe Ansatz method [35]. The eigenvector for a system of 
N particles can be parametrized by N complex numbers Zk with 



N 



(x<)= J2 M{zi})U z ?(k) ( 22 ) 



ct£Sn k — 1 



where the sum is performed over all the permutations a of N elements and x< = (xi, . . . , Xjy) is the vector containing 
the ordered positions 1 < xi < X2 < ■ ■ ■ < < L oi the N particles on the ring. The coefficients Ap and the N 
numbers are fixed by imposing the eigenvalue condition W^Vi^ = ^ 4>\ , up to a global normalization constant 

( s) 

of the vector 6) ' . 



The N Bethe equations resulting from the Ansatz (22 1 are given for all j £ {1, ... , N} by: 

4 = (-i)^ 1 TT p e ' + g e "'^-(p + gfr (23) 

kyy peS + qe SzkZ i ~ (P + 1) Zk 

with the eigenvalue: 



N 

pe 



fc=i 



+ qe s z k - (p + q) ) . (24) 

z k 



One can read 1 18 L 123] for easy derivations of these equations and eigenvalues, as well as for the expressions of the 
coefficients A a ({zk})- 

This set of non-linear equations has many solutions but is expected to contain all the eigenvalues and thus the 
maximal one Ai(s). For s = 0, the matrix is stochastic and the first eigenvalue is known to be Ai(0) = 0; 

moreover, the corresponding eigenvector has constant coordinates since all the configurations are shown to have the 
same probability. In terms of the Bethe roots, it corresponds to the Bethe roots Zk going all to the same value 1 as 
s 0. A perturbative study of Ai(s) around s = for all sizes L done in jSSJ HU [M| has already given the expression 
of the cumulants of the current for the ASEP for various asymmetry regimes and the first finite-size corrections. In 
particular, the attempt of resummation of the large deviation function f(j) from the cumulants has shown that a 
singularity (finite convergence radius) occurs for the same critical values j c and s c as in the hydrodynamic approach. 
Nevertheless, the perturbative expansion does not allow to understand what happens in the jammed phase, nor at the 
phase transition. 

The translation invariance of implies that the eigenvectors (jif^ are eigenstates of the translation operator. It is 



indeed the case as it can be seen from the construction (22 1 and the eigenvalue is the product Z\ . . .Zjsr- One can check 
from (23) that one has (z\ . . . zn) l = 1, which corresponds to the periodic boundary conditions. In the case of the first 
eigenvalue Ai(s), the eigenvector <fi^ is non-degenerate and all its entries are positive (Perron- Frobenius) ; moreover 
it has to be an eigenvalue of the one-site translation operator with eigenvalue exp(2i7rn/L), n £ {0, . . . , L — 1}: the 
reality and the positivity of the entries of 0j select together the integer n = for Ai(s) and one has to satisfy 

N 

TI z k = 1 (25) 
ft=l 

for the first eigenvalue Ai(s). The Gallavotti-Cohen symmetry [57] corresponds to the change of variable e s h-> (q/p)e~ 3 
and Zk H> 1/Zfc. 
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III. BETHE EQUATIONS FOR THE WASEP 
A. Polynomial representation 



The form (23 1 of the Bethe equations is not suitable to a simple analysis of the weak asymmetry regime s ~ j/L 
and p/q ~ 1 + 2v/L. The change of variable 



Uk\/qjp + 1 



inserted in (23) leads to the following set of equations on the u/c's: 



../ {1 Vj, e s 



n 



puj - qu k 
quj - pu k ' 



(26) 



(27) 



The change of variable is such that Uj goes to as s — > and z k — > 1 simultaneously and is defined as long as 
Zfc 7^ e s {p/q) and u k { — \/p/q, —\Jqjp\. The Gallavotti-Cohen symmetry corresponds to the change u k i-> l/u k . 
One notices that the interaction terms in the r.h.s. of (23) and (27 1 are now rational fractions with numerator and 
denominator of degree 1. The polynomial Qn(U), of degree N, is defined as 



N 



Qn{U) = IJ(tf-Ufc), 



(28) 



k=l 



so that the set of N equations ( 27 1 is equivalent to the divisibility of the polynomial 

L 



e sL p N 



U 



1 



-u 



p 



1 



Q 



N 



-u 



by Qn(U), since they are both vanishing on the Uj, that are supposed to be distinct in the generic case. Thus, there 
exists a polynomial Rl,n(U) of degree L such that : 



R L ,N{U)Q N {U) = e sL p N [UJZ + I] Q N 



-V 



P 



i Q 



N 



-u 



<1 



(29) 



The r.h.s. of this equation is linear in Qjq but is non-local in U since it implies the computation of Qn at other points 
than U. There does not exist, to our knowledge, any systematic method to solve algebraic equations such as (29 1 for 
all values of L and N. Only asymptotic solutions for L, N — ¥ oo are available under some additional hypotheses on the 
asymptotics of the Bethe roots, as well as some relations with Askey- Wilson polynomials [33]. 
The translation condition (251 can be written in terms of Qn(U) for the eigenvalue Ai(s): 



n ( p\ n QN{-\Jq/p) 



97 Qn(-Vp/i) 



= 1. 



Eq. ( 29 ) for U = gives the normalization condition 

R L ,N{Q) = e sL p N + q N 

of the polynomial Rl,n- 



(30) 



(31) 



B. The resolvent and its properties 

The purpose of this sequence is to translate the algebraic divisibility property of Qm for finite N and L into analytic 
properties when L, N — ¥ oo with N/L — > po and the weak asymmetry scaling. Eq. (29 1 can be rewritten : 

L 



Rl,n{U) 



2{Jpq) N (\ + Uy 



-sL/2 



1 



sL/2+N In y/p/q 



1 



1 



p/qU ] QN{{g/p)U) 



+ U I Q N (U) 



i p -sL/2-N In yJvTq 

l + VqfpU\ Q N ((p/q)U) 
\ 1 + U J Qn(U) 



(32) 
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The ratios Qn(V)/Qn(U) can be rewritten with the so-called resolvent Wl,n(U): 

Qn(V) 



Q N (U) 6XP ( L Su WL > N(u)dU )> (33a) 
1 N 1 

W^ (t0=B ^ (33b) 
fe=i K 

The function Wl.aK^) has JV poles located at the itfc's, is holomorphic everywhere else on C and decreases as po/U as 
\U\ —> oo where N/L = po. 

The resolvent technique presented here is standard in the study of the thermodynamic limit of Bethe equations and 
a simple overview of this technique in a similar algebraic context is presented in [22 . 

In most systems solved by coordinate Bethe Ansatz, the Bethe roots that correspond to the ground state Ai 
condensate on a finite number of smooth curves or degenerate points in the large size limit, although this property 
has not been shown rigorously. Numerics for the WASEP are in agreement with this observation |2(il |3~4"] . This 
condensation property translated for the resolvent function Wl,n(U) says that it converges to a limit W(U) that has 
a finite number of cuts or poles in the complex plane, which correspond to the smooth curves or accumulation points 
formed by the itfc's, and W(U) is holomorphic everywhere else. 

It leads us to the following assumption, on which all the computations of the next sections rely: 

Assumption III.l. In the large N,L limit with N/L — > p , Ls — > 7 and L(p/q — 1) — > 2v, the resolvent Wl,n(U) 
associated to the maximal eigenvalue converges to a well-defined limit W(U) with a finite number of poles and cuts and 

holomorphic everywhere else. In the same limit, the function [\+u) L conver 9 es t° a well-defined limit r{U) for 

all U 7^ — 1, that is holomorphic on C — { — 1} and may have a singularity at U = —1 that can be an essential one. 



If this is true, then (32) has a limit given by 



r{U) = cosh ( I + pv + v ^-jj - 2vUW(U) V (34) 

In the weak asymmetry limit and under the previous assumption, the three equations (30), (24 1 and (31 1 become 

7P + 2p^ + 2//T^(-l) = 0, (35a) 
LA 1 ( 7 /£) i ^>A 1 ( 7 ) = 2is 2 (w(-l)-W'(-l)) 1 (35b) 

r(0) = cosh (Z+pv)- ( 35c ) 
By construction, the resolvent W(U) takes the generic form 

"■poles (i) "cuts r (j)/ \ 

i=l 2=1 JF 3 

where p( % ' > is the fraction of roots condensed on the pole and p^\z) is the density of roots condensed at the 
point z of the oriented cut Tj. The resolvent W(U) can be related to the Stieltjes transform of the spectral measure 
often introduced in random matrix theory, where the eigenvalues are replaced by the Bethe roots in the present case. 
One shall notice that the existence of poles for W(U) is compatible with the assumption that the Uk's are distinct for 
finite L and N: finite size effects may split immediately the roots around the poles U^K The density p^\z) and the 
fractions p^ can be extracted from the knowledge of the resolvent W(U) through: 

= lim 01- U®)W(U), (37) 

U-H7W 



\U) = lim f W(U - ie) - W(U + iej), VU eT,-. (38) 

e->0+,e along Tj ZITT \ I 

In order to make the notations easier to read, it is useful to define for every U on an oriented cut Tk the two operators: 

A(->/(t/) - lim (f(U-ie)-f(U + ie)) (39a) 

e— >0+,e along T fc V / 

A' + V(P) = lim (f(U-ie) + f(U + ie)) (39b) 

£— >0+,£ along F k \ / 

and one has A^f(U) = 2inp^(U) for U £ Tj. 
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C. Characterization of the poles and the cuts 



As a limit of rational fractions with only one multiple pole at U = — 1, the function r(U) which appears in (34 1 



is analytic on C — {— 1}. This implies that the r.h.s. of ( 34 1 can not have any cuts nor poles, excepted at U = — 1 



From the change of variable between the z^'s and the u^'s in the limit L — > oo, it follows that the u^.'s can not be 
accumulated at U = — 1 and thus W(U) is well-defined at U = — 1 (it can also be seen on ( |35a[)). The only possible 
pole of W compatible with the analyticity properties of r(U) is at U = 0. In the decomposition (|36|, at most one pole 
remains at U — with a residue p > 0. 



Equation ( 34 1 can be rewritten as 

r(U) = cosh$(Z7) (40) 

where $(J7) contains the same amount of information as W(U) in its cuts and poles. This function $(£/) has a simple 
pole at U = — 1 with residue —v and has a finite limit for U — !• oo. Applying the operator A^~) defined in p9| to both 
sides of the previous equation gives for U on a cut of $(?/): 

= sinh (A ( - ) $([/)/2) sinh (A (+) $(£/)/2) . (41) 

This means that on any point U of a cut of one has either the condition 

A (-) $(Z7) = 1mm, meZ, (42) 

or the second condition 

A (+) $([/) = 2i7rn, n G Z. (43) 

The two functions A^ ± ^$(/7) are continuous along one cut, so the value of the integers morn have to be constant along 
a cut, until one can jump from one condition to the other if there is a point where they are satisfied simultaneously. 
This remark is the key point of the identification of the phase transition in section |IV B| 



Conditions of the first type (42 1 give directly through (38) the density pV'(z) on the cut of W{U) where it is 
satisfied: 

P«HU) = ^L. (44) 



The root density on the cuts for which (43 1 is satisfied is more difficult to obtain; it can be performed by solving 
Cauchy principal value equations as in |2l| but this tool will be useless in the simple case described in the next section. 

Finding the set of cuts with the integers rrij and corresponding to the largest eigenvalue Ai(s) is not a simple 
task. The next sections present the choice that gives the Gaussian part of the large deviation function and a candidate 
for the jammed phase in agreement with numerical simulations [261 134] . 

IV. LARGE DEVIATION FUNCTION OF THE CURRENT 
A. Gaussian part of the current fluctuations 

It is known from the stochasticity of the matrix that the rescaled eigenvalue Ai(7) has to vanish when 7 — > 
and it corresponds to all the Bethe roots Zj going to 1, i.e. all the u/s going to 0. The resolvent is then given by 
W(U) = p/U when 7 = 0, which has a single pole at U — and no cuts. 

As soon as 7 ^ 0, cuts have to be opened around some points in the complex plane. Numerics [26; '31] tend to show 
that the largest eigenvalue A 1(7) corresponds to the choice of a single cut T for a discontinuity ( |43| with n = 0. The 
study of this cut is presented below and one gets easily convinced by similar computations that any opening of a 



second cut of any type (42 1 or (43 1 is impossible within the analyticity constraints on cf> and r, once assumption III. 1 
is accepted. 

Let us introduce the extreme points ao and a\ of the cut T, the polynomial P2(U) = k(U — ao)(U — ai) such that 
P%{— 1) = 1 and the function 

fa(U) = \fPW) (45) 
which has the same cut T as &(U) going from a and a\ and the asymptotics 4>2(U) ~ U as U — > 00 in order to define 



completely the sign of the square root. Then, if one assumes (43 1 with n — on the only cut T, the product $(U)4>2(U) 
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does not have any cut anymore on the complex plane by construction: the only remaining pole is at U = — 1 and the 
asymptotic behaviour is &(U)(p2(U) ~ 0(U) at +00. One knows from complex analysis |3S] that it must take the form 
S 2 (U)/(U + 1) where S2 is a polynomial of degree at most 2. Moreover, r(U) has no other pole except at —1: $([/) 
cannot diverge for U — > ao,i and Pa(C^) has to divide S2{U). One gets finally: 



$([/) 



U+l 



(46) 



Only three parameters {n, clq and a{) are left. The properties (351 of W can be translated in terms of the function 
$(£/): 



$([/) 

$([/) 
cosh$(0) 



U + l 



A. 



Ai(7) 



(U +l)+(D((U + 1) 2 ), 



U— >oo 



= cosh 



with the definitions of the two coefficients 



A-i 



7(1 - 2p) 



1 
2 



1/(1 



and the values of the three unknown quantities k, ao and ai can be computed easily. One obtains 

P 2 (U) = v 2 - 2vA_ x (U + l)+A 2 ao (U+ l) 2 . 



(47a) 

(47b) 
(47c) 

(47d) 
(47e) 

(48) 



One also checks that (47c) is simultaneously satisfied. Computing the eigenvalue Ai(7) is then an easy task, through 
an expansion of &(U) at the next order around U = —1: 

Ai(7) 



A 2 - A 2 



p(l - p) 

2 



2-yv 



(49) 



It corresponds to the result expected from the macroscopic hydrodynamic theory (JJ. 

Although it has not been necessary for the computation of the eigenvalue, the shape of the cut of the Bethe roots 
can be extracted from <&([/)• The resolvent W(U) is given by: 



W(U) 



2vU 



U 



pv 



-$([/) 



The point U — belongs to the cut if and only if 



Po 




U + l 

v < 7 < —2pv. For 7 > —2pi>, the residue at is given by : 

if 7/1/ > 0, 
p if -2pv < 7 < 0. 



(50) 



(51) 



The shape of the cut depends on the sign of 7, as shown on figure [2] For 7 > 0, the polynomial Pi^U) has two complex 
conjugated zeros and the cut T can be obtained by considering A^W(U) = pr(u)e ie ^: the angle 9(u) gives the 
local direction of the cut and pr(u) gives the density of the roots. For — 2pv < 7 < 0, the zeros and the cut are real 
and the density along the path is simply given by : 



Pr(u) 



2u{u + l)' 



(52) 



For —v < 7 < —2pu, the pole U — belongs to the cut and the density contains an additional Cauchy principal value 
(the pole at U — creates a local hole in the density, which produces this principal value) . It is interesting to note that 
the eigenvalue, as well as the function <&, is not sensitive to the change of behaviour of W(U) around but depends 
only on the polynomial PzQJ) whose coefficients remain continuous as 7 varies. 

The shape of the cut and the density of the roots along the cut gives a curious parallel with random matrices : the 
GUE ensemble of matrices with Gaussian entries has a spectrum given by the Wigner semi-circular law. Here, after 
change of variable, the Bethe roots have a semi-circular law whereas the fluctuations of the densities are expected 
(from the hydrodynamical limit) to be Gaussian. No understanding of this potential parallel is known however. 
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Figure 2: Condensation of the Bethe roots for different values of 7 in the Gaussian phase for p = 1/2: for y/u > 0, the cut is 
symmetric under z — > z whereas it is real for ■y/u < 0. This second case 7/1/ < also contains a pole at 0. The densities along 
the cut in the Gaussian and the travelling wave phase are illustrated in figure [3] 



Pr(x) 



Pr(a 




x i-> l/{2vx) 



ao cii 
Gaussian phase 




x h-> l/(2vx) 



a b 61 
Travelling wave 



Figure 3: Density of roots along the cuts in the Gaussian and the travelling wave phases (solid lines) : compared to the first 
case (section IV A I, the cut in the second case (section [IV C \ is split in three parts and involves elliptic integrals. 



B. The phase transition 



As 7 varies, one checks that the creation of a small new cut of type (42 1 or (43 1 somewhere else in the complex plane 
is impossible within the analyticity constraints on $ and r. However, there exists another scenario of modification of 



the cut obtained for the Gaussian part that gives a larger eigenvalue than (20 1 and (49 1 for 7 beyond a critical j c . 

The previous cut was obtained by considering the case (43 1, which sets the value of A^ + )$(C/) to on the cut. The 
computation of A^~^$(C/) on the same cut T gives: 



A (_) $([/) = 2: 



U + l 



(53) 



This quantity is equal to 2mm for any m £ Z only on a discrete set of points on the cut. The trivial case m = 
corresponds to the edges of the cut where the density vanishes. The case m = ±1 is satisfied only if v^J-P 2 {U)-tt{U+1) 
vanishes on the cut, which correspond to a local density satisfying p cu t(U) — 1/{2vU). One checks that this is possible 
only if 



Pi 1 - P) t-2 



7T~ 

y 



(54) 



This threshold is precisely the one obtained from the linear stability analysis in the macroscopic hydrodynamic theory 
[7] but it is derived here only from complex analysis in a non-perturbative way. 

At the transition point p{\ — p)(7 2 + 2^v) — —ir 2 , there exists a single point b on the cut on which both constraints 
( |42| and (43 1 are satisfied. For p(X — p){j 2 + 2^v) < —n 2 , the double point b splits into two points bo and b\, between 
which one can satisfy either (42) or ( [43] ) (see figure [3]). The second choice corresponds to the Gaussian case studied 
previously, whereas the first choice leads to another eigenvalue, which is computed in the next section. 

The phase transition occurs necessarily when 7^0 and all the cumulants of the current are not affected by the 
phase transition at the leading order in the system size. The phase transition can be seen in the cumulants only 
through the divergence of the finite-size correction, as described in [211 124| . 
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C. Beyond the phase transition: a single hybrid cut 



Beyond the transition, the candidate cut for $(J7) is made of four singular points: the two extreme points ao and a\ 
of the cuts and the two intermediate points 60 and 61 where the cut condition switches from (43 1 to (42 1. The cut T of 
$ can be split into three parts : Tq (from ao to 60) and I\ (from b\ to ai), on which (43 1 is satisfied with n = 0, and 
r_ between them (from 60 to 61), on which (|42| is satisfied with m = 1 (see figure (3|. 



As a generalization of the Gaussian case described previously, one can introduce the polynomial Pa(U) = (U 
ao)(U — bo)(U — bi)(U — ai) and the associate square root function <f>4,(U) defined by 



(55) 



and whose first (resp. second) cut joins ao and bo (resp. b\ and ai) and coincides with the part T (resp. Ti) of the 
cut of $([/) with type (43). The arbitrary sign of the square root is fixed by imposing 4>i(U) ~ U 2 as U — > 00. 
The derivative $'([/) satisfies the following cut conditions: 



A(+)$'([7) = on T and Ti, 
A(") $'([/) = onL. 



(56) 



The function $'([/)</> 4 (J7) is thus holomorphic on C — { — 1}, has a pole of order 2 at U = — 1 and is 0(1) as U — >• 00 
by construction of W(Z7): it has to be equal to T2 ([/)/([/ + l) 2 where T2 is a polynomial of degree 2. Around —1, we 
have the expansion $'([/) = v/(U + l) 2 — A/i/ + o(l) and thus $'([/) is given by: 



> 4 (-i) + (g + + 4>'U-W + 1) 2 /2 - + 1) ; 

(u + i) 2 Mu) 



(57) 



where A is the corresponding eigenvalue. The function $(£/) can be written either as an integral of $'({/) or by the 
following formula obtained by analysing the cuts and poles of ^(U)/4>4(U): 



u + i 



61 



2 + 1 1 



dz, 



(58) 



which involves only the polynomial P4 and satisfies directly both cut conditions ( 42 1 and ( 43 ) . Matching the previous 
expression (58 1 with the required asymptotics (47 1 leads to: 



bi 



bo 



bo 



Z + l 

<j) A {z)' 



dz, 



7 (1 - 2 P ) 



bi 



+ vp 
Ai(7) 



4 (-l) + ^(-l)(z + l) 

(2 + l)0 4 (^) 

i 4 (0)(z + l) 



z4>a{z) 



dz, 



rbi 
Jbo 



(z + l) 2 </) 4 (z) 



dz. 



(59a) 
(59b) 
(59c) 
(59d) 
(59e) 



There are four independent constraints (59) for four unknown quantities ao, bo, b\, a\ and the eigenvalue Ai(7) is 
completely determined. The sign chosen for (59dl (left free by the condition 47c I is set in order to be in agreement 
with the density of roots condensed in zero in the Gaussian phase (51 ). 



The expression of W(U) can be recovered from ( |58[ ) and the definition of <3> (C/) . It gives the density of Bethe roots 
along the cut: 



1 

2vu 



p T {u) = < 



4tti>u(u + 1 



61 



z + l 



60 



for u € [6 ,6i], 



dz for 14 g [ao,6o] U [61, a{\, 



(60) 



12 



A(7) 



15 



10 



-10 



-15 



Atrayelling (7) 

AGauss(7) 
-7T 2 /2 






-25 



-20 



-15 



-10 



Figure 4: Legendre transform of the large deviation function for p = 1/2 and v — 10. This plot shows that the hybrid cut with 
both A' + ' and A' - ' conditions gives a larger eigenvalue as soon as A(7) < — 7r 2 /2. 



and the density of roots po condensed at u = 0: 

7 + 2pi> 



Po = 



Av 



1 

27 



60 



VFf£(o)T* + i 



dz. 



(61) 



The case corresponding to a cut including and a negative fraction po can be treated as in the Gaussian case by 
considering a Cauchy principal value of the quantities written above; it does not lead to any additional phase transition 
nor any relevant non-analyticity. 



D. Partial correspondence with the hydrodynamic description 



Numerics show that the expression of the eigenvalue we obtained coincides with the one obtained from the 
hydrodynamical approach described in appendix [A] However, the exact mapping between the polynomial P4. which 
satisfies the set of equations (591, and the polynomial R4, which satisfies the set (A5), seems to be non-trivial (at least 



more complicated than just an homographic transformation) . 

For the particular case p = 1/2 however, the mapping is exact and simple. For a half-filled system, the speed v of 
the travelling wave in appendix [A] is v = and the profile is symmetric under g — > 1 — g (particle-hole duality) : we 
obtain the values Ci = and g\ = \ — gi (with the notation of appendix |A|. This particularity emerges as well in the 
Bethe Ansatz approach and corresponds to an additional symmetry of the polynomial P^U). 

In the Gaussian case, the polynomial P 2 is given by 



P 2 (U) = - V 2 U + 



7 + v 
2 



[U + lf 



(62) 



and is invariant under the additional symmetry P2(U) — » U 2 P2(1/U). 

In the non-Gaussian case, the polynomial Pa{U) shares the same additional symmetry P 4 {U) U i P 4 {l/U) = P 4 {U). 
The equation for A_i fixes one of the coefficients such that Pi{U) can be written as: 



P A (U) = P U 2 +p'U(U + l) 2 + (U + l) 4 = (U+ 1) 4 P 



u 



(63) 
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The two coefficients p and p' satisfy 



7 + v 



1/60 



1 



bo y/Pi{u) 

1/6(1 v^HO 

bo V^lM 



(in, 



(64) 
(65) 



which are equivalent to (A5dl with v = and (A5a) after the change of variable g = U/(U + 1) and a different 
normalization of P 4 ; the eigenvalue is the same in both cases. 

The full correspondence between both approaches for p 7^ 1/2 is however not clear and requires further investigation, 
as well on the correspondence between the polynomials P4 and R4 as on the relation between the real travelling wave 
profile and the curves of the Bethe root. 

A second interesting explicit correspondence corresponds to the limit v — > 00 with the scaling 7 = 1/7; this limit is 
expected to coincide with the TASEP large deviation function as already explained in [7 . The divergences that appear 
in equations (59 1 imply that do f fro an d a\ 4 b\ in this scalings. In terms of the characteristics of the cut, it means 
that the two parts with A^ + ^$ = shrink to zero and only remains the cut with = 2i-K. Two methods lead to 

the same result: either one extracts carefully the asymptotics of the all right terms (59) or one constructs directly <E> 
with the prescribed limit behaviour. One obtains for any U not to close to the cut 



$([/) - v 



a 



Qj-i 
17 + 1 



In 



U-bi 

U-bn 



•0(1/) 



(66) 



Satisfying the renormalized conditions (47) leads to a_i = 0, ao = 72 + (1 — p) and 

'1 + V 



In 



l + br 



-jp, In 



"7 



It leads to the scaling of the eigenvalue: 



lira 



Ai( 7 ) (l- e w)(l-e^( 1 -")) 



^-+00,7/^-4-7 V 



(67) 



for all 7 < 0, and one recovers the asymptotics of the TASEP as explained in [7J. It is interesting to see that the 
coexistence of high and low density domains in the TASEP behaviour corresponds to A( _ ) cuts, whereas the gaussian 
fluctuations around a flat profile correspond to A^ + ) cuts and it would be of interest to understand this correspondence 
from the eigenvector ( 22 1 itself. 



E. Existence of other phase transitions : a short discussion 



A question left open in jJJ is the existence and the nature of other phase transitions in the WASEP that may 
overcome the phase transition described in the previous sections. 

The standard scenario in Bethe Ansatz (as well as in random matrices) is a condensation of the Bethe roots on 
curves or isolated points in the thermodynamic limit N, L — > 00 with N/ L — > p. These condensation curves correspond 
to cuts of the resolvent W. Phase transitions occur when the number or the topology of the curves vary: splitting of a 
curve into two, emergence of a second curve, closure of a cut, etc. 

The scenario presented here is different: no change occurs in the number of cuts nor in their topology (and the 
pole of W{U) at is irrelevant) but a phase transition still occurs due to the change of type of the Riemann surface 
associated to the function $>(U). Such a phenomenon does not appear in random matrices (GUE with a potential for 
example) since the resolvent of the eigenvalues satisfies a quadratic equation much simpler than the cosh equation 
(40 1. The Riemann surface of the argcosh functions gives further possibilities, which are well described by the set of 



conditions ( 42 ) and ( 43 1 
As alreac 



y announced in section IV A it is not possible to open other cuts of type ( 43 1 . One remarks first that 



the function $'([/) has cuts only of type (43 1. If we suppose the existence of n cuts with end points a,j and a' and 



construct the polynomial p2 n (U) whose zeros are these end points, then f(U) = $' ' (U) y/p2n(U) has no cut if one makes 
the cuts of the square root coincide with the ones of 4>'(?7). Thus, / is analytical on C except at U = — 1 where it has 
a pole of order 2. Moreover, f(U) = 0(U~ 2 ) for U —> 00. The only choice corresponds to f(U) = R n (U)/(U + l) 2 
where R n (U) is a polynomial of degree n. The analyticity properties of $ set then R n (U) = except for n = 1, which 
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corresponds to the Gaussian case. Opening a cut with (42 1 only is also impossible. Thus, there appears to be no phase 
transition due to a change in the number of cuts. 

However, within the one-cut assumption, one could imagine a hybrid cut as in the previous section but with the 



condition (42 1 for \m\ > 2 or even a hybrid cut made of successive paths alternating between the two conditions (42 1 



and (43 1, generalizing the scheme presented in the previous section. Checking if it is possible or not requires a precise 



study of the integrals of the type ( 59 1 that would emerge with a polynomial of larger degree: it is beyond the target of 
the present paper but it would be very interesting to check it rigorously. 

The discussion of the existence of several curves presented in this section is valid only for part of the eigenvalues 



(more precisely, the ones for which the assumption III. 1 is valid). When the asymmetry p — q is reduced from the 



PASEP to the WASEP, the largest eigenvalue as well as some others follows the scaling (20) for which assumption III. 1 
holds for the Bethe roots. The other part of the spectrum does not have to satisfy this assumption and this scaling 
form: with another well-adapted scaling, the Bethe roots may still condensate on some curves in arbitrary number. 

V. CONCLUSION 

The exclusion process in the weak asymmetry regime p — q = 0(1/L) shows a phase transition in the large deviation 
function of the current that was first predicted in [7]: this paper presents an alternative analysis of this transition with 
the Bethe Ansatz toolbox. 

Different types of results have been obtained in the present paper and each of them leads to some new questions. 
From the point of view of integrable systems, the first interesting point is the shape of the condensation curve of the 
Bethe root in the thermodynamic limit considered in this paper: to our knowledge, it is the first example of a single 
cut with a piecewise-defined continuous density. Its origin is the non-trivial cosh equation satisfied by the resolvent 
instead of the standard quadratic equations. 

One of the interest of the scaling p — q = 0(1/ L), which is unusual for spin chains, is the existence of an hydrodynamic 
limit that leads to partial differential equations that are classically integrable whereas the microscopic Markov transition 



matrix satisfies quantum integrability (see section II B \. The present example gives another example of transition from 
quantum to classical integrability. 



The main progress that leads to the solution presented here is the formulation ( 32|33| of the Bethe equations that is 



suitable to the scaling limit p — q = 0(1/L). Variations on this procedure should also lead to interesting results |36] 
for other weak asymmetry regimes such as the KPZ scaling regime p — q = 0(1/V L), which is studied actively for 
different geometries |26l l2"81-30j . Higher order perturbative expansions should give access to the finite size corrections 
of the large deviation function. It would be particularly interesting to characterize more precisely the phase transition 
(exponents, scaling of the finite size corrections, etc.). 

Some open questions also emerge from the present paper. Although the same integrals that determine the large 
deviation function of the current appear in both the Bethe Ansatz and the hydrodynamic approaches, the precise 
correspondence between the contour of the Bethe roots and the density profile of the travelling wave is still missing: it 
would be interesting to understand the relation between them and to establish a "dictionary" between corresponding 
notions in both approaches. A corollary would be the correspondence between the two fourth-degree polynomial P4 
and i?4 that emerge in the two approaches. 

The last new information that comes with the Bethe Ansatz is the knowledge of the eigenvectors, which also describe 
in principle the spatial correlations. Computing correlation functions from the Bethe Ansatz is known to be a difficult 
problem in general; nevertheless, the macroscopic hydrodynamic description tends to show that the macroscopic 
correlations should have a simple description. It is probable that the coordinate Bethe Ansatz form of the eigenvectors 
should have a simpler limit than usual in the present weak asymmetry regime but it has not been possible to study it 
rigorously yet. 
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Appendix A: Optimal travelling wave profile of the variational principle 



The normalized Legendre transform (|9| of the large deviation function of the current is obtained from the variational 
principle ( 17 1: 



A1X7) = sup 7 J — lim 



J,p(x,r) 



Z[0,T\ (j, P) 



T->c 



T 



(Al) 



where the fluctuating profile p and current j satisfy the conditions (15) and (16) with jo — J. The optimal travelling 

(A2) 



wave studied in [7] corresponds to the choice of a particular type of time-dependent profiles 

p(x,t) =g(x-vt). 



Ai(7) = sup ( 7J- 

J,g,v 



One obtains j(x, t) = J + v(g(x — vt) — p) from ( |15[ ) and (16 1, and ( Al I becomes time-independent: 

1 UJ + v{g(x)-p)-u<j{g{x))f g'jxf 1 \ 
MM) 8a( g (x))\ J' 

The previous functional is invariant by translation of the profile g and its optimization yields (see [7]): 
g\xf = 4(J + v(g(x) - p) - vcj{g{x))f - 4(d + C 2 g(x))a(g(x)) = 4R ( /' v \g(x)) 



(A3) 



(A4) 



where C\ and C2 are two constants that are fixed once the extremal densities g\ and gi of g(x,t) over [0, 1] are fixed 
by the two conditions: 



dx 



g{x)dx = 



-dg, 



4 J ^(g) 

92 n 



» ^j4 Lv \g) 



dg. 



(A5a) 
(A5b) 



The change of variable x v-+ g[x) for the two intervals over which g is monotonic yields: 



Ax (7) = sup ( 7 J - hc x + C 2 p) ~ 

J,v \ ^ 

The optimal values of J and v are then obtained implicitly by: 



y-2 



1 



7 + 1/ = j 



= j 



1 



« a(g)y/R{ J ' v \g) 



dg + u 



.9 - P 



« cr{g)>jR?' v) (g) 



z dg + w 



« a{g)^t%) 

(9-p) 2 
« cr(g)jR?' v) (g) 



dg, 

dg. 



(A5c) 

(A5d) 
(A5e) 



For p = 1/2, the speed is w = and the correspondence with the Bethe Ansatz calculation presented in section IV D 
is given by J = v^fpi and Ci — 2^J = p%. 
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